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We study magnetic domain walls in double exchange ma- 
terials. The domain wall width is proportional to the square 
root of the stiffness. In the double exchange model the stiff- 
ness has two terms: the kinetic energy and the Hartree term. 
The kinetic energy term comes from the decrease of the tun- 
neling amplitude in the domain wall region. The Hartree term 
appears only in double exchange materials and it comes from 
the connection between band-width and magnetization. We 
also calculate the low-field magnetoresistance associated with 
the existence of magnetic domains. We find a magnetoresis- 
tance of 1 — 2%. The magnetoresistance can be considerably 
larger in magnetically constrained nanocontacts. 

PACS number 71.10.-w, 75.10.-b, 72.10 

Mixed valence compounds of the form Lai-xAxMnO^ 
{A being Ca, Sr or Ba) have recently been shown 
to present extremely large (colossal) magnetoresistance. 
10,^ In these materials strong Hund's interaction between 
the charge carriers and the manganese ions leads to a 
strong coupling between the electrical resistivity and the 
magnetic state. For 0.1 < x < 0.5 and low temperatures, 
the system is metallic and presents ferromagnetic order. 
As the temperature increases the system becomes insula- 
tor and paramagnetic. The magnetic transition occurs at 
a x-dependent critical temperature Tc ~ SOOif. Colos- 
sal magnetoresistance occurs near Tc and in presence of 
magnetic fields of several Teslas. 

In the Lai-xAxMnO^ compounds, the electronically 
active orbitals are the Mn d orbitals, and the mean d oc- 
cupancy is 4— X. A strong ferromagnetic Hund's rule cou- 
pling aligns all electron spins in the Mn d orbitals. The 
Mn ions form a simple cubic lattice of lattice parameter 
a. The cubic crystal symmetry splits the d orbitals into 
a t2g triplet and an Cg doublet. Three electrons fill up 
the t2g levels forming a core spin of magnitude S — 3/2 
and the rest of the electrons, 1 — x per Mn, go to the Cg 
orbitals. 

Ferromagnetism in these materials is explained by the 
Double Exchange (DE) mechanism, in which the elec- 
trons get mobility between the Mn ions using the oxy- 
gen as an intermediate. This conduction process 
is proportional to the electron transfer integral and due 
to the strong ferromagnetic Hund's rule coupling it is 
maximum when the two cores spins involved in the pro- 
cess are parallel and it is zero when they are antiparallel. 
Because the alignment of spins favors electronic motion, 
the ferromagnetic ground state maximizes the electron 
kinetic energy. When the temperature increases, the DE 
model undergoes a phase transition towards a paramag- 



netic state. In this phase the core spins are randomly 
oriented and the electron kinetic energy is minimized. In 
the paramagnetic phase these materials behave as elec- 
trical insulators. |^ 

Large low-field magnetoresistance has been observed in 
ferromagnetic Lai-^AxMnOs compounds with different 
structural discontinuities. |^-^ These effects are associ- 
ated with a lack of oxygen at the interface which produces 
an antiferromagnetic ordering at the interface and breaks 
the DE mechanism. |]l0| 

Below Tc these materials may contain magnetic do- 
mains separated by domain walls (DW's). Domain walls 
produce a resistance to the electrical current, and for un- 
derstanding low-field magnetoresistive effects in DE met- 
als it is important to know the width and the resistance of 
DW's. Domain wall magnetoresistance effects have been 
recently observed in itinerant ferromagnet systems such 
as Co, Ni and Fe [ pd]JT2| and also in colossal magnetore- 
sistance perovsquites. [|13[ In reference [13] it is obtained 
that in Lai-xCaxMnO^ the resistance of a domain wall 
is 8 X 10~^'^Um'^, a quantity that the authors argue is 4 
orders of magnitude larger that one might expect based 
in DE models. Recent theoretical works have studied the 
ballistic and diffusive transport through domain walls in 
itinerant ferromagnets, |Q however for DE systems only 
the transmission in one-dimensional models have been 
studied. (l5). 

In this paper we study magnetic domain walls in DE 
systems. In the DW the direction of the Mn spin changes 
from to TT over a region of width Lw- In this region 
the Mn core spins are misaligned and the tunneling am- 
plitude between Mn ions along the DW is reduced. This 
loss of kinetic energy in the DE model is the equivalent 
to the loss of exchange energy in the Heisenberg model. 
The chemical potential in the system is fixed by the mag- 
netic domains which have all the same hole concentration 
X. The reduction of the bandwidth in the DW region 
with respect of the surround magnetic domains produces 
a change in the density of electrons in the DW region. 
This effect cost a lot of Hartree energy and the system 
prefers to create dipoles at the edge of the DW's. In this 
way the Mn ion levels change and the local charge is not 
modified in the DW. The shift of the energy level of the 
Mn ions modifies the cost of creating a DW and there- 
fore its width. This effect is new and occurs due to the 
DE mechanism. In the first part of the paper we charac- 
terized this effect and study how it affects the width of 
the DW. In the second part of the paper we study the 
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ballistic transport through DW's in DE systems. 

Microscopic Hamiltonian. We are interested in a hole 
concentration in the range 0.1 < a; < 0.4. For this doping 
range the Cg orbitals are degenerated and for simplicity 
we consider that the Hamiltonian is degenerated in the Cg 
orbital index. Also in our model the Mn spins are treated 
as classical. For temperatures below Tc and in the limit of 
infinite Hund's coupling, the electronic properties of the 
Mn oxides are described by the nearest neighbor tight 
binding Hamiltonian, 
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Here C^^ creates an electron at site i, in the orbital a 
and with spin parallel to the core spin at site i, and the 
hopping amplitude is given by mm 
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where 9i and (pi are the angles which characterize the 
orientation of the core spin at site i. 

Long-wavelength functional. In a perfect system all 
core spins are parallel and the electron kinetic energy gets 
its maximum value. Now consider a modulation of the 
core spin direction characterized by a vector M; = Si/ S. 
In the long-wavelength limit the modulation can be de- 
scribed by a continuum unitary vector field M(r). This 
modulation produces a decrease of the value of the hop- 
ping amplitude, and therefore a kinetic energy loss. For 
smooth modulations the local loss of kinetic energy is pro- 
portional to (VM)^ and the difference in kinetic energy 
(KE) between the uniform system and the modulated 
system is given by 
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where p^^ is the KE stiffness of the system. This term 
is the equivalent to the term describing the exchange en- 
ergy loss in the Heisenberg model. Note, however that 
in the DE system p^^ is associated with the loss of elec- 
tron kinetic energy due to the spatial variations of M. By 
diagonalizing the Hamiltonian Eq. 1 with different bound- 
ary conditions we have calculated the value of p^^ for 
different hole concentration in the system. In Fig.l we 
plot p^^ as a function of x. Note that we have a degen- 
eration 2 associated with the Cg orbitals and we measure 
the hole concentration with respect half-filling. When the 
electron concentration is zero (a; = 1) the band is com- 
pletely empty there is not KE in the system and p^^ = 0. 
When the band is half filled (x = 0) the KE is maximum 
and p^^ gets its maximum value. 

Besides the KE, there is also a Hartree (H) term which 
contributes to the long-wavelength functional. The elec- 
tron density, n ,in different regions of the sample is fixed 



by the chemical potential p which is obtained from the 
value of n in the regions of constant magnetization (mag- 
netic domains). The tunneling amplitude only depends 
on the relative orientation of the Mn spins and therefore 
n and p get the same values within all the magnetic do- 
mains. These regions represent basically all the system 
and they can be considered as reservoirs for the elec- 
trons. The total charge in the system should be zero and 
a background of positive charge equal to 1 — x is assumed 
to exist in the sample. 

In regions with modulated magnetization, VM ^ 0, 
the tunneling amplitude is reduced and the bandwidth is 
narrowed. Since p is fixed by the reservoirs, the decrease 
of the bandwidth would produce a change in n in these 
regions. This breaks local charge neutrality and it would 
cost a lot of Hartree energy. The system prefers to create 
dipoles in order to shift the Mn ion energy levels in such 
a way that local charge neutrality is recovered. The en- 
ergy shift is negative (positive) for electron concentration 
smaller (bigger) than half-filling. For smooth modulation 
of the magnetization, we have calculated from the Hamil- 
tonian Eq.(l) the energy shift required for keeping local 
charge neutrality. We find that the shift is proportional 
to (VM)^ and therefore the local charge neutrality con- 
strain gives a contribution to the energy of the system of 
the form, 
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In Fig.l we plot the value of /3 as a function of X, as 
obtained by solving the Hamiltonian Eq.(l) with different 
boundary conditions. In principle we should add also a 
term describing the energy cost of creating the dipoles 
responsible for the energy shifts, however we find that 
this energy is much smaller than AE^ . 

Finally we need to know the energy term corresponding 
to the constrain which creates the magnetic domains. In 
general this term has the form. 



AE^ = J d^rf{6) 



(5) 



For uniaxial crystals f{9) = _R'sin^(6') being K the 
anisotropy constant. For domain walls created by mag- 
netic constrictions [p^|jl3| f{0)=gpBHS/a^ cos6'sgn(x), 
being H the magnetic field. 

Adding the different contributions we get the func- 
tional 



F^^J d^r{VMf + J dhf{e) 
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being p = p^^ + 213. 

Domain wall width. A DW along the x direction has 
the general form M — (0, sin 6*, cos 6*) and cos9{x) is ob- 
tained by minimizing F with the boundary conditions. 
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M ±z at X —I- ±00. For uniaxial crystals the opti- 
mum form of the DW is 



cos = — tanh 
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with — \J p/'2'K. For a magnetic constrain the form 
of the domain wall is given by 
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and L\Y = "^y PjJI^TTs' "'"^ both cases Lw represents the 
width of the DW, and it is proportional to the square root 
of the total stiffness p. The effect of the Hartree term to 
the DW width depends on the electron concentration: 
for electron concentrations smaller than half filling f3 is 
negative and the Hartree term prefers thin DW's. On 
the contrary for electron concentrations bigger than half 
filling P is positive and the Hartree term likes wide DW's. 

For a given value of Lw equations (7) and (8) have 
essentially the same form and for simplicity we use ex- 
pression (7) for describing a domain wall. 

Transport through a domain wall. Now we calculate the 
ballistic conductance, G, associated with a DW. From the 
difference between the conductances of a perfect system 
and a system with a DW we can evaluate the low field 
magnetoresistance associated with the alignment of the 
magnetic domains. A DW modulates the magnetization 
only along x direction and the Hamiltonian is invariant in 
the y — z plane. Therefore the conductance of the system 
can be written as 



G{p)= / dE'gME')n^''{t^~E') 
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here n [E) is the two-dimensional density of states per 
unit area at energy E, and gmiE) is the conductance 
of a one dimensional system at energy E. Expression 
(9) can be interpreted as the sum of the conductance of 
all the one dimensional channels. The one dimensional 
conductance between the sites 1 and N is written as, 



(10) 



being D{p) = At^ — p/ (27ri^) the local density of states 
at an edge of the isolated lead, and GN^i{fi) the Green 
function connecting the sites 1 and N of an infinite one 
dimensional system. The factor 2 in Eq.(lO) corre- 

sponds to the Cg degeneracy. The DW is fully contained 
between the sites 1 and N. 

In order to calculate the effect of the DW on the trans- 
port properties it is necessary to know the effect that the 
magnetization, M = (0, sin 6*, cos 0) has on the electron 
Hamiltonian, There are two effects: the modulation of 
the hopping amplitude along the x direction and the shift 



of the Mn ion levels. Both effects are obtained from the 
form of the DW, Eq.(7). 

From Eq.(9) we evaluate the magnetoresistance. 
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Here Gdw and Go are respectively the conductance in 
presence and in absence of the DW. MR represents the 
low-field ballistic magnetoresistance associated with the 
presence of magnetic domains. In Fig. 2 we plot, for dif- 
ferent hole concentrations, MR as a function of Lw- For 
small values of Lw, MR gets rather large values (> 10%). 
For Lw > 20a, the magnetoresistance is always smaller 
than 1%. In the inset of Fig. 2 we plot the the one- 
dimensional conductance, guj as a function of the energy 
for a DW with Lw=5a. The asymmetry with respect to 
zero energy, is due to the Hartree energy level shift. We 
see that giu is suppressed mainly at the band edges. This 
is the reason why the MR is bigger for small concentra- 
tion of electrons, where the Fermi energy is close to the 
band edge. 

Discussion. We now discuss the application of our 
results to manganese oxides. First we want to know 
the width of the DW produced by crystal anisotropy. 
The width is determined by the stiffness, p, and by 
the anisotropy constant K. K can be obtained ex- 
perimentally from neutron scattering p9[-[2l|], by mi- 
crowave absorption and by studies of the resistance 
saturation with magnetic fields. Q From these experi- 
ments we estimate an anisotropy constant in the range 
2 X 10^ - 5 X 10^ for x = 0.3 this implies a DW 

width in the range 10a < Lw < 30a and a ballistic MR 
of 1 - 2%. 

In the case of DW's created by magnetic constric- 
tions ||l^ , Lw depends on the applied magnetic fields. 
For X = 0.3, we obtain Lw — 26a/ VH being H the 
external magnetic field in Teslas. This corresponds to 
a rather large DW. Therefore we expect that the DW 
width should be determined by the crystal anisotropy 
and ~ 1 — 2% magnetoresistances are expected. This 
value is similar to the obtained in reference [13], assum- 
ing 10a < Lw < 30a. However we do not know what 
is the contribution of diffusive processes to MR. Within 
the simplest Born approximation, diffusive effects are de- 
scribed by the reduction of the bandwidth, and for the 
values of Lw considered this produces also a 1% magne- 
toresistance. 

From the inset of figure 2, we can say that the mag- 
netoresistive effects should be much more important in 
geometrically constrained domain walls. Recently mag- 
netoresistances of 200% have been observed in nanocon- 
tacts of itinerant ferromagnets. |l^,^. We expect simi- 
lar effects to occur in DE materials confined to reduced 
dimensions. 

Summary. We have calculated the width of domain 
walls in double exchange materials. The width is propor- 



3 



tional to the square root of the stiffness. The stiffness has 
two terms; the kinetic energy and the Hartree term. The 
kinetic energy is the equivalent to the exchange energy in 
the Heisenberg model and it comes from the decrease of 
the tunneling amplitude in the domain wall region. The 
Hartree term appears only in double exchange materials 
and it comes from the connection between band-width 
and magnetization. We have also calculated the low- 
field magnetoresistance associated with the existence of 
magnetic domains. We have found a magnetoresistance 
of 1 — 2%. The magnetoresistance can be considerably 
larger in magnetically constrained nanocontacts. 
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FIG. 1. Variation of a) p^^ and b) /3 as a function of x, 
as obtained from Hamiltonian Eq.l. Note that x is measured 
with respect half filling. 
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FIG. 2. Variation of MR as a function of Lw for differ- 
ent hole concentrations. The inset shows the one-dimensional 
conductance as a function of the energy, for a DW with 
Lw=5. 
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